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Xl 



LIST OF SYMBOLS (continued) 

a EMG standard deviation 

pressure differential in tube of length I, N/m 

^ degrees of angular position measured counterclockwise about 

the optical axis (Figures 2. 11 and Z. IZ) 

T] kinematic coefficient of viscosity, m /sec 
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SECTION 1 
INTRODUCTION 



This report presents the development of a new technique lor 
enhancing the fidplity of flight simulators. The name giv^ to the tech- 
nique and to the system which applies the technique is BIOCONAID, an 
acronym for Bionic Control of Acceleration Induced Dimming. 

When an aircraft is subjected to a change in its velocity vector an 
acceleration is experienced by the occupants. The cyDmmon situation is 
that of a pilot changing the direction of motion of the aircraft through a 
steep turn or pull-up. Then, if the maneuver is coordinated, the pilot 
experiences an acceleration terding to puil him tighter against his seat 
as if he weighed more. This acceleration also affects the blood pres- 
sure distribution in his body. The blood pressure su^Dplying the eyes and 
the brain is reduced and the pilot may experience gray-out or black- 
out if the turn is continued or the pull-up becomes steeper. 

In a combat situation, the pilot n^edi to continue the turn or pull-up 
and retain his vision* He may perform straining maneuvers (described in 
Appendix C) to increase the eye-level blood pressure and fight the black- 
out tendency. Such straining represents expenditure of energy which 
leads to fatigue. The "G-suit" is a protective garment which exerts pres- 
sure on the lower extremities, thereby restricting the pooling of blood to 
the lower extremities during high-acceleration maneuvers. Pilots of 
high performance aircraft wear this garment to increase their tolerance 
to acceleration. 

Modern flight simulators, which are widely used for combat crew 
training, provide visual scene dimming as the siriialated G-load increases. 
Also provided are G-suits which realistically react to the G-load. Missing 
however, is the pilot's need to perform the beneficial straining maneuver 
with its attendant energy expenditure. Also missing is the characteristic 
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loss of vision progressing from the peripheral field inward. The 
BIOCONAID System provides the means to simulate both of these missing 
elements. 

The BIOCONAID System consists of four integrated software models 
plus hardware which senses and processes pilot straining information for 
input to the computer. The computer software models are: 

• Cardiovascular Model (predicts effect of G-level on eye-level 
blood pressure) 

• G-suit Model (predicts beneficial effect of G-suit at the input 
G-level) 

• Straining Model (predicts the beneficial effe-.t of straining) 

• Visual Field Model (predicts the visual scene brightness and 
the limits of the visual field under acceleration) 

The complete system also includes an electromyographic (EMG) subsystem 
(signal processor). 

S<»ction 2 of this report includes a discussion of fl^e development of 
each of the models mentioned above plus the EMG Subsystem (EMG strain- 
ing signal). Also treated is the integration of the entire system. Sections 
3 and 4 cover development and documentation of the algorithms and stra'n- 
ing experiments. 

Section 5 gives conclusions and recommendations. Brieily, the 
BIOCONAID System performs as planned. A usable, smoothed EMG signal 
is produced by pilot straining. Either the abdomen (external oblique) or 
the back (latissimus dorsi) muscles are suitable as electrode sites. 
Neither extraneous motions nor electronic noise seems to affect the system 
adversely. AH software algorithms perform satisfactorily. Full implemen- 
tation of the Visual Field Model was not accomplished during the contract 
period due to equipment limitations. 

Because possible detrimental cardiovascular effects of sustained 
straining in a 1 -G environment may exist, it is recommended that some 
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medical research be conducted to identify system limitations or restrictions. 
It is also recommended that the BIOCONAID' System, with possible limitations 
or restrictions, be evaluated in an Air Force simulator with experienced 
pilots. 

Appendix A gives a test plan for system evaluation which might be 
accomplished in the F-4 simulator at Luke Air Force Base. Appendix B 
gives a complete listing of the computer programs developed during this 
effort. Appendix C gives a discussion of the M~l/L-1 straining maneuvers 
and Appendix D describes the special Tustin solution method used for the 
cardiovascular model. 
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The human organism is a highlv complex system of interconnecting and 
interdependent systems. Its responses, in general, must be considered non- 
linear and extremely difficult to model. However, by restricting the number 
and levels of variables in a specific situation, it is often possible to model 
complex systems successfully. 

The purpose here is real-time modelling of the human response to G- 
loading* To simplify the problem, it is assumed that only z-directed (normal) 
acceleration will be considered and that the only variables influencing accelera- 
tion are G-levei, G-suit protection, and pilot straining. The complex human 
system can then be partitioned into four tractable models: the cardiovascular 
model, the protective garment model, the straining model and the visual scene 
model. The relationship between these models is shown in Figure 1. The 
result of each of these blocks is a contribution to the pseudo blood pressure 
variable called ^^protection variable" (PV). In all cases, the contribution to 
this blood pressure variable is either derived directly from the acceleration 
profile or is a sc;condi.,ry result of that profile and the operation of life support 
equipment. Note that a set-point value (one-G blood pressure at eye level) is 
summed with the PV values (see Figure 1). A normally distributed random 
variation, labeled "n" is also included to account for variations with pilots on 
different days. 

The cardiovascu ar response model is represented by a linear transfer 
function which corresponds to the acceleration response of the pilot's blood 
preasure. One prime goveming factor in pilot response to acceleration is the 
onset of grayout and blackout. These visual problems are directly related to 
the blood pressure at eye level. This model output provides a dynamically 
correct signal which is equivalent to the nominal loss in eye level blood 
pressure for an unprotected human undergoing any acceleration profile, G(t). 
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The straining simulation mode! accounts for the G tolerance enhance- 
ment which is afforded by a properly executed M-) or L-1 straining maneu- 
ver. The purpose of these straining maneuvers is to increase the eye level 
blood pressure in order to maintain vision. Proper performance of the 
maneuvers requires that the abaominal and upper torso muscles be tensed 
isometrically and that expirations should be made against a closed glottis. 
The result is an increased intrathoracic pressure and increased blood pres- 
sure at the eye. Proper application of the straining maneuvers also results 
in ihe appearance of myoelectric signals on the skin surface. These low 
voltage electric signals are directly associated with muscle contraction 
(Bigland and lippold, 1954; Brody, Scott, and Balasubir amanian, 1974). 
These biologically derived signals are then processed by the model to gener- 
ate a straining contribution to the protection variable, PV , which represents 
the increased blood pressure due to the M- 1 or L-1 maneuver. 

The pressurized G sui^J^s *a protective garment used to increase thf 
individual's tolerance to +G^. The suit uses pressurized bladders to press 
against the legs and lower abdomen. The external pressure inhibits pooling 
of the blood volume in the lower extr emities* thus insuring a better supply- 
to the heart during acceleration. The suit must be inflated l^y the G valve 
to a predetermined level for the suit to be effective. The simulation model 

accounts for the required pressure level and uses the actual suit to provide 

a . 

the necessary dynamics. The suit pressure is compared with the normal 
flight pressure schedule and a protection variable contribution is generated 
by the model. 

The dynamic visual field model was developed to be readily implemented 
in a simulation system. The visual field model uses the eye level blood pres- 
sure predicted from tne simulated G profile and the other protection variable 
contributions as inputs to produce dynamically responsive signals which pre- 
dict the expected visual field of a pilot undergoing that G profile, 

The system presented on the following pages represents the results of 
partitioning this complex physiological system into models. Each 
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model is explored and developed in detail based on a common protection 
variablepP\4 which is related physiologically to system blood pressure. 
The blood pressure models are then finally combined in a systematic 
paradigm which provides the driving values for the visual field response 
model. 

2, 1 THE CARDIOVASCULAR RESPONSE MODEL 

Human tolerance to long term acceleration is normally measured 

in terms of visual loss (blackout) and unconsciousness. Both of these 
tolerance end points are related to the ability of the cardiovascular system 
to deliver oxygenated blood at adequate pressure to the retinal and cerebral 
regions. The acceleration causes a changing blood pressure profile in 
the human such that the effective pressure vertically abcve heart level is 
decreased and the pressure below heart level is increased. There is 
therefore a lower perfusion pressure at* eye level. The distribution of the 
blood in the body also changes as the acceleration pools blood in the lower 
parts of the body and lungs. There is therefore less available blood to 
circulate and a lower oxygen content because the longb do not operate as 
efficiently. 

There are two cardiovascular systems which are dynamically involved 
in the process of blood pressure maintenance while the human is under- 
going ' +G acceleration. The hydrostatic system is responsible for the 
reduced retinal perfusion pressure at the eye and eventual loss of pressure 
at the cerebral level. The orthostatic system is related to blood pooling in 
the lower body with a concomitant reduction of venous return to the heart. 

The cardiovascular system has feedback systems which are affected 
by blood volume and pressure. The feedback systems attempt to regulate, 
the pressure and flow characteristics of the cardiovascular system. One 
of the primary pressure sensors for this system is located in the carotid 
artery. Any change in pressure at the sensors results in a regulatory 
operation which attempts to return the pressure to a preset value. For 



example, a pressure drop at the carotid sinus initiates a neurologic reflex 
action which changes the heart rate and output so as to increase the overall 
pressure leveL This section develops the p^tionale for a cardiovascular 
transfer function model which describes the acceleration response of the 
blood pressure delivery s/stem. 

A radically new model for the cardiovascular system wcs not the 
intent of this system. Rather a reasonably accurate transfer function 
based on physiological and physical principles as well as experimental 
observation was sought. The major requirements were that the transfer 
function model retain the major response characteristics of the cardio- 
vascular system, be implementable in a real time digital simulation system, 
and retain factors which are identifiable with measurable physiological 
factors. 

A number of models of the cardiovascular system have been developed 
(Gillingham, Freeman, and McNee, 1972; Knapp, Randall, Evans and 
Marquis, 1978; Koushanpour and Spickler, 1975; and Miller and Green, 
1973). The model developed by Gillingham (Gillingham, ct al. , 1972) is 
used in the present effort. This linear model uses a transfer function of 
the form: 

K fl + as) J'V (s) 

^^^^ , 2 G(s) 

1 + b^s + b^^s 

Values selected for a^, b^ and b^ are compromise values from the 

literature and Gillingham's response curves* The break frequency for 

the lead compensator factor, (1 * a^s), was selected at 30 mHz (millihertz). 

2 

The denominator factor (1 + b^s + b^s ) was selected to give complex roots 
with a natural frequency of 70 mHz and a damping ratio of 0.7. These 
characteristics result in time responses and Bode plots which agree 
reasonably well with experimental measurements. To obtain these 
characteristics we find that a^ = 5. 31, ^bi 3. 23 and = 5. 17. 




In order to evaluate the system gain constant, K^, Rogers, 1978, 
treated the artery from the aortic arch to the eye as a column of blood 
under che hydrostatic action of the G-loading. This colunm is offset from 
the Z-axLS by an angular displacement of 9 degrees due to the angular 
inclination of the pilot's seat. The system gain determined in this way 
- -21.4, compares very closely with the gain derived from experi- 
mental blackout results (Burton, Leverett and Michaelson, 1974), When 
the seat inclination is included, becomes -Z1.4 cos The final 
transfer function is shown in Figure 2. (Note that G(t) here is dynamic 

G- loading - G -1.) The solution technique used for generating PV 
z c 

from G(t) is described in Appendix D. 



Figure 2. Block Diagram of Cardiovascular System Response 

2.2 THE PROTECTIVE GARMENT MODEL (G-SUIT MODEL) 

The history of G- suits dates from World War II and has been well 
reviewed in recent monographs (Burton, et al., 1974; Eraser, 1966; 
Gaur and Zuidema, 1961; and Howard, 1965). The wrap-around C$lJ-3r/P 
cutaway-type of anti-G suit, presently used by the U.S. Air Force, 
improves blackout tolerance in the +G direction by about 2G above 
resting tolerance (Leverett, Whitney and Zuidema, 1961). However, the 
mechanism by which this beneficial effect takes place is not entirely 
understood. Investigations have shown that the protection afforded by 
the suit is not associated with any increase in cardiac output (Lindberg 
and Wood, 1963) nor with any increase in eye-level blood pressure 
(McCally, 1970). The G-suit attempts to prevent pooling of the blood in 
the lower extremities of the body during plus-G maneuvers (Leverett, 



G(s) 



-21. 4 cos 9 (1 i 5. 31s ) 

2 

1 + 3. 23s i 5. 17s 
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et al., 1961). Thus, it is not possible to model the physiological system 
dynamics pertaining to the G-suit until more knowledge is gained about the 
system. 

The approach taken in the present research ffort is to use an actual 
valve plus a G-suif- on the pilot to p^^ovide ^he suit dynamics and to model 
only the known beneficial effect-- the ZG improvement in blackout 
tolerance with a properly operating suit. The valve will be activated by 
the computer-generated G signals. 

The normal suit pressure versus G-lcading schedule is shown in 
Figure 3» If the actual suit pressure, at a given G-level, falls in the 
cross-hatched area labeled "Acceptable Region", then partial protection 
is being provided. Full ZG protection tis afforded when the suit pressure 
reaches *'he right-hand edge of the "Acceptable Region". At higher suit 
pressures no additional protection is provided and at suit pressures to the 
left of the "Acceptable Region" no_ protection is provided. 

A linear variation of protection within the acceptable region is 
assumed which has a width of approximately 80 mm Hg. The simulation 
procedure is to obtain the difference tween the actual suit pressure and 
the stcindard curve pressure {right hand edge of acceptable region, 

Pigure 3). This difference is labeled AP^. Then the protective garment 

G 

protection value contribution is generated as follows; 

if tP^ <-80 then PV^ = 0 (2) 
G G 

if-80i:AP^<0 then " ^ ^ ^ ^G ^ ^ 

if 0 S AP^ "then PV_, = 4Z. 8 mm Hg (4) 

G G 

Note that 4Z.8 mm of Hg is equivalent to 2G protection since the cardio- 
vascular system gain is -Z1.4. The simulation mod^l is presented in 
block diagram form in Figure 4. 
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Figure 4. G-Suit Model for Generation of Protection Value PV 
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2. 3 THE STRAINING MODEL 

Human tolerance of acceleration may be improved by the use of pro- 
tective measures and devices. Straininij; and gruntinj in tight turns was 
practiced by German pilots prior to World War II to improve their 
tolerance. This technique gradually evolved into the M-1 m^neu\er. At 
the same time it was discovered thart bending forward also impro\ ed the 
tolerance of pilots in high acceleration maneuvers. The effects of posture, 
restraint, and body position on G tolerance have been well worked out and 
have been reviewed in detail (Davson, 1969; Fraser, 1966; and Gaur and 
Zuidema, 1961). 

Since the early days of flying, the value of muscular straining as 
protection against the effects of +G acceleration has been recognized. 
Pilots observed that blackout or grayout could be postponed if they grunted 
or screamed and tensed the skeletal musculature during acceleration. The 
M-1 maneuver, defined as muscular straining with expiration against a 
partially closed glottis (see Appendix C) is a most effective voluntary pro- 
tection against the circulatory effects of +G acceleration (Lindberg, Sutterer 

z 

Sutterer, Marshall, Headley, and Wood, I960; White, 1958). The physi- 
ological basis for this protection has been ascribed to its effect on 
increasing the arterial pressure at eye level during acceleration. As 
forced expiration ij instituted the resulting increase in intrathoracic 
pressure is transmitted directly to the aorta. There is an immediate 
increase in peripheral arterial pressure which is only slightly less than 
the increase in intrathoracic pressure. This maneuver affords up to 2.4 G 
protection against blackout (Burton, et al., 1974). The circulatory 
mechanisms described were identified by Rushmer in 1946. The "pr* ^sure 
raising" approach has theoretical difficulties in that blood pressure is 
controlled within rather narrow limits by a number of effective servo- 
mechanisms, the baroreceptor system being the bedt known. As a matter 
of course the elevated intrathoracic pressure can result in a decreased 
venous return, lowered cardiac output and therefore a lower blood 
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pressure at the eye* The carotid sinus reacts to the initially increased 
blood pressure and also attempts regulation. Both of these effects are 
nmininnized by a cyclic application the M- 1 or L~ 1 procedure with a 
repetition period of 4-5 seconds. 

The skeletal muscle is controlled by signals wh^ch are transmitted 
to selected motor ^mits of a muscle through the motor neurons. The force 
generated by the muscle is the result cl both the freqaency of firing of 
motor units and the number of motor units which are recruited. The 
muscle tension is accompanied by electrical signals (EMG) which can be 
detected by suitable electrodes on the skin surface. The electrical signal 
exhibits the characteristics of its primary source in that it represents a 
weighted sum of motor unit activations (Shannon, 19V5). There is therefore 
a correspondence between the EMG and muscle force (Bigland and 
Hppold 1954). 

The EMG signal can be considered to be the summation of a large 
number of individual potential waveforms each resulting from the firing of 
a motor unit. Thus we can assume that the EMG signal is Gaussian with 
a mean of zero. If such a signal is the input to a full wave rectifier, it has 
been shown (Davenport and Root, 1958) that the mean output is proportional 
to the standard de viation of the input. Furthermore, the standard deviation 
has been shown (Shannon, 1975) to be the appropriate measure of the EMG 
signal to use for control purposes. This reasoning leads to selection of 
the rectified, smoothed EMG signal for determining a straining protection 
value . 

The simulator model requires the presence of two signals to generate 
the straining protection value, PVg . The two signal3 represent the prime 
factors which are present in a properly executed straining maneuver. The 
EMG signal is generated by the straining subject and processed to provide 
in intermediate protection value PVj^^. The protection afforded by muscle 
straining is then modified by the effects of the breathing maneuver (M-l)» 
The protection factor ^c^nerated by the breathing pattern is critically 
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controlled by the straining and breathing rate. An improperly perforired 
M-1 can actually reduce protection and result in a detrimental pressure 
drop. 

We assume the thoracic cavity to be represented by a flexible 
cylinder surrounded by a musv^le girdle. The muscle straining signal 
representing muscle force around the thoracic volume is then assumed to 
- bear a linear i^elation to the internal pressure increase at low straining 
levels. At higher levels a maximum pressure is achieved and increases 
in muscle straining are no longer effective. The pressure straining 
relationship used in the model is shown in Figure 5 as a curve with a 
saturation level related to the maximum pressure rise in thu thoracic 
cavity. 

The effectiveness of the straining maneuver is also governed by the 
repetition rate. The straining protection value PVj^ is therefore modified 
by a function dependent upon maneuver repetition period ^^^^R^* Experi- 
mentally, K (t ) is maximum at 3-5 seconds according to Gillingham 
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(Appendix C). Shorter time periods do not allow sufficient time for the 



responses due to baroreceptor feedback and reduced venous return to the 
heart. The function representing the variation of with maneuver 

repetition rates is shown in Figure 6. The straining simulation block 
diagram is shown in Figure 7. 

2.4 THE VISUAL FIELD MODEL 

Rogers, 1978 discusses the structure and functioning of the eye as it 
forms the basis for a usable model of the visual field under acceleration. 
The effect that acceleration has on the visual apparatus is observed in 
terms of tunnel vision, greyout and blackout. During the periods of 
grayout there are a^so decreases in visual acuity and brightness contrast 
detection ability. Tunnel vision means that less of visual perception 
occurs first in the peripheral areas of the visual field and progresses 
jnward. Although there are man/ factors related to the anatomy, psychology 
and physiology of the human which are responsible for these effects, the 
structure of the eye is a primary factor, and it provides a suitable model 
for the effects of acceleration on the visual field. 

The primary cause of the effects of accelerationoon vision is believed 
to be the reduction of blood supply to the rods, cones and cells in the retinal 
layers. ITiis has been confirmef^. by various investigators " '^wis and 
Duane» 1956; Newsom» Leverett and Kirkland, 1968). The v s supplied 
by the opthalmic artery which has two branches in thf^ cye itself. One 
branch, the retinal artery, enters the eye trough the optic nerve and pro- 
vides th*^ primary source of blood to the retinal layers. The other branch, 
the ciliary artery, supplies the outer layers of the eye and the ciliary 
processes. Figure 8 shows the retinal artery entering through the optic 
disc and branching out to supply the entire retina. 

The structure at the retina is idealized as a neiwork of cylindrical 
hydraulic pipes with a central supply at the entrance of the optic nerve. 
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Figure 6. Effect of Repetition Period on Straining Protection Values 
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Figure 7. Straining Model for Generation of Protection Value (P^^) 
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Figure 8. Topography of the Eyeball 

from W. J. White in Reference Zl^ 
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Figure 9. If a steady laminar flow is assumed through a cylindrical tube 
the flow, Q, will vary with the 4th power, of the radius of the tube and the 
pressure differential AP. Also the flow will vary inversely as the length £ 
of the tube and the kinematic viscosity of the fluid T|, This is stated in 
Poiseuille's law as: 

Q.ILlilAP) ' (5) 

The blood pressure drop across the cardiovascular system is mostly 
across the arterioles (Berne and Levy, 1972). At higher levels of pressure 
the retinal arteries supply the arterioles with sufficient pressure that the 
entire retina is perfused with a continuous flow of freshly oxygenated blood. 
As the pressure at the retinal artery drops, the flow must decrease in 
accordance with Poiseuille^s Law, and less blood, therefore less oxygen, 
is available to replenish the retina. 

The effective pressure (P^) in the recinal artery is the difference 

R 

between supply pressure P^^^ and the interocular pressure Pjq pressure 
differential of 10-20 mm Hg is required to maintain eyeball rigidity). 

P = P - P 
R ae lO 

As an initial approximation, the retinal supply network is assumed 
to be a network with the supply pressure distributed linearly across the 
network. That is at any point along the supply. 

dP 

r 

; — = constant 
dx 

where x is the distance in degrees along the arteriolar bed, P is the 

r 

corresponding differential pressure along x,and the slope is a function of 
the supply pressure at the optic disc. 

The fam^ily of curves describing the pressure gradient across the 
retina is a series ol straight lines with slope 
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where P„ represents the current retinal arterial pressure at the disc and 

is the outer circulation limit of the retina corresponding with the ora 
seratta. For any value of P^^^ the retinal pressure as a function of degrees 
from the disc is given by: 

In a small blood vessel ^ pressure differentia) of about 15-20 mnn Hg 
across the vessel wall is required to maintain an open flow channel (Berne 
and Levy, 1972). When this transmural pressure is not maintained, the vessel 
collapses, and the blood supply is cut off to all points beyond. This critical 
pressure differential is now labeled P 

CRIT 

The threshold value for any x (called x ) in iTigure 10 occurs when 

Lj 

P falls below P . Solving for x in terms of P at the P 

T i^Kii jLt ae CxvL r 

^CRIT=-(H7^)''L^<^ae-^I0' 

/P - P - P \ 

„ / ae 10 CRIT \ 

-L-S\ P - P^^ j (8> 

^ ae 10 ^ 

For a nominal value of = 110^ representing the outer circulation limit, 

x^ is given by: (let P^^ = 20 mm Hg, P =17 mm Hg) 

L ^ ' 10 ^ CRIT ^ 

110(P - 37) 

^L^ P ^^20 
ae 

The quantity represents arc length measured from the center ot 
the optic disk. The optic disk is located on a visual field map approximately 
10^ to the right of the optical axis and 10^ below the horizontal right eye 
and to the left and below for the left e>e. The peripheral edge of the model 
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Figure 10. Poiseuille's Law Applied to the Retina 
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is determined by a polar projection of the retinal field. The coordinate 
origin of the polar map is taken as the visual axis and the optic disc is 
taken as the flow center. Defining il as degrees of angular location 
measured counterclockwise about the optical axis (Figure 11), the field 
contribution of the right eye under acceleration is then given by the solution 
of: 

2 2 2 

r - 2r r cos (d + 4^) + r , x (10) 
d a L 

where r (degrees) is the sensed visual limit at ^ » offset distance 

(degrees) of the disc from the visual axis and d is the angular offset 
(degrees) of the disc from the horizontal axis. A similar expression is 
found for the )eft eye. 

The visual field model can now be extended to a binocular field. 
Figure 12 is a binocular field map showing the coincident foveal areas as 
the visual center, an^ the optic disc for each eye in a manner that depicts 
the observed field from inside out. The outer lines depict the outer edges 
of the peripheral vision. Lowered blood pressure supply in each of the 
modeled retinas causes the field to collapse as concentric circles with 
centers at each of the optic discs as shown by the dashed lines. Thus the 
visi^al field "ollapses toward a somewhat oval shape with the vertical field 
having a smaller visual angle than the horizontal field. This shape is in 
agreement with Gillingham's experimental results (Gillingham and 
McNaughton, 1972). 

At some point after the visual field has collapsed^ the central vision 
is lost. iilingham, et al., 1972, shows that central vision remains after 
the peripheral field has collapsed. The operating model is therefore 
modified to allow for central vision by incorporating a central visual cone 
which is responsive to the eye elevel protection value. Mechanization of 
the model utilizes a binary decision to determine presence of the central 
field. The decision point is chosen at a nominal 10 mm Hg below that of 
the peripheral field critical point. 
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In the actual mechanization of the visual field model, the following 
data are used: 

visual cone for foveal vision - 5^ (Gillingham, et al., 1972) 

r^ - 10^ (for both left and right eyes) 
- . 175 radians 

d - 10^ (right eye) measured clockwise from the right horizontal 
visual axis 
= . 175 radians 

d = 170^ (left eye) measured as for right eye 
= 2. 967 radians 

110 (P - 37) 
ae 

x^ ^ — 20 degrees (Equation 9) 

ae 

1.920 (P - 37)/P - 20) radlans 
ae ae 

P in this equation becomes PV, the generated protection 
ae 

value, which is equivalent to eye-level blood pressure under 
acceleration. 

The above information is used in Equation 10 which is repeated 

below: 

2 2 2 

r - 2r,r cos (d + ^) ^ r^ = x_ (10) 
d a ij 

All angles in degrees are converted to angles in radians and the equation 
is solved for "r", which is the visual limit at angle t about the visual axis. 

In either field: 



r^r cos(d + *)±/r_^^ cos^ (d + * ) - ( - x^) (ID 
d d a L, 

In the right field; 

r ^ . 175 cos (.175 + ^) 
r 



✓ ? P V - 3 7 2 

.031 cos (.175 + t) - .031 + 3.686 ( > '^^^ 
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In the left field: 

r = . 175 cos (2. 967 + ^) 
I 

/2 PV-3 7 2 

.031 cos (2. 967 + '];) - . 031 + 3. 686 ( ) (13) 

Thus, using Equations 12 and 13, the sensed visual limit may be 
calculated for any angle ^If and for any protection va^ue PV. These equations 
are applied only to the point where all peripheral vision has collapsed and 
only the central (foveal) vision remains. The visual cone for foveal vision 
is selected as 5^ (Gillingham, et al., 1972). In the software, the collapse 
point is reached at a value of PV which leaves the quantity under the radical 
in Equation 12 positive: 

. ..031 . 3. 686(^^)^ = 0 ^'^^ 
PV-20 

which gives, 

jr^V ^ 38. 72 (15) 

The relationship between retinal blood oerfusion, and the ability for 
contrast detection in a visual scene is apparently an oxygen transpo^-t 
phenomenon. The observed visual dimming aiid field coUap-^^ model could 
possibly be extended to include these effects. Chambers and Hitchcock, 1963, 
have documented the relationship between contrast detection threshold and 
the G field, and White (1958), has documented the effect of lowered visual 
acuity due to reduced oxygen content of blood in the retinal artery. Using 
the CV model for the response of eye level blood pressure to acceleration, 
the contrast and acuity results could be added to the visual dimming and 
field collapse modeU 

In the present simplified mechanization, it has been assumed that the 
brightness and contrast, at any given point in the visual field, are linearly 
related to PV, the eye level blood pressure as influenced by acceleration. 
When PV is known, the Visual Response Algorithm produces, in addition to the 
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visual field limit, signals appropriately scaled to be used for controlling 
both brightness and contrast at any point in the visual field. 

The algorithm presently includes control of a center and two side 
visual TV screens with separate brightness and contrast signals. Collapse 
of the peripheral vision will res^ilt in the two side screens being dark 
while collapse of the foveal vision will result in all three screens being 
dark. Prior to collapse, the screens would exhibit a linear dimming. 

2.5 THE EMG SUBSYSTEM 

Figure 1 shows the EMG subsystem between the pilot and the strain- 
ing model. The task of the EMG subsystem is to sense the raw EMG signal 
from the pilot's skin, amplify it, process it; and present it to the digital 
computer. 

The expanded block diagram of the EMG subsystem is shown in 
Figure 13. Three electrodes (one is a reference ground) pick up the 
l>ilot*s EMG signals. The bio-amplifiers add the two electrode voltages 
ind provide a gain of 100. A hi-pass filter elmiinates low frequency drift. 
An^plification up to 1000 then readies the signal for full-wave rectification. 
Finally, a two- stage low- pass filter provides a smooth usable signal for 
use in the simulator computer. 
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Figure 13. EMG Subsystem Block Diagram 

Except for the electrodes and the bio-amplifiers, the EMG sub- 
system, as shown in Figure 13, was patched on an EAI TR-20 Analog 
Computer. The patchboard wiring diagram is shown in Figure 14. 
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2. 6 SYSTEM INTEGRATION 

The required dynamic model of pilot visual response to G can now be 

z 

assembled as a summation of the separate models. The completed system 
algorithms provide the means for implementation m aircraft training simula- 
tors. The integrated system is shown in Figure 15. 

Pilot commands enter the simulator aircraft dynamics causing 

changes in the aircraft velocity vector. The resulting g-loading G(t), is used 
to drive the cardiovascular model • In addition the G(t) signal causes the 
pilots G-suit ^) to inflate according to a preselected schedule. The effect 
of the proper suit inflation is to increase the value of the protection variable 
within certain Limits. The sum of these signals ^) is added with a computer 
generated random signal related to individual variations and with a set value 
representing the eye-level, nonaccelerating, blood pressure (120 mm Hg). 
The resulting PV signal drives the visual field model. The instrument panel 
and the view screen are then dimmed according to the output of the model ^) . 
The pilot senses the size and brightness changes and has the option of reducing 
his aircraft maneuver intensity (g-loading) or increasing his g-tolerance level 
by performing a straining maneuver (M-1 or L-1). If he initiates such a man- 
euver, the EMG subsystem processes the myoelectric activity from his skin 
surface^ providing a smoothed EMG signal ^) to the straining model. The 



straining model then generates a protection value, ' which is depen- 

dent on the straining level and on the straining interval. This PV is added 

s 

with the other effects and thus serves to increase PV and to increase the 
visual field limits and brightness. 
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Circled numbers refer to Figure 15. 
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Figure 15. System Integration Diagram 



•SECTION 3 

ALGORITHM DEVELOPMENT AND DOCUMENTATION 



3. 1 INTRODUCTION 

Two FORTRAN algorithms have been developed from the physiolog- 
ical models presented in Section Z. These algorithms are: 

1. The Dynamic Visual Response Algorithm 

2. The EMC Algorithm 

The Dynamic Visual Response Algorithm consists of a FORTRAN 
coded implementation of the Cardiovascular Response Model, the Protective 
Garment Model, and the Visual Response Model. The EMG Algorithm 
consists of a FORTRAN coded implementation of the EMG Straining Model. 
A description of each algorithm is presented along with its flow chart. 
In addition, a program listing combining the two algorithm^ and the program 
results with several input profile are presented. 

3.2 DYNAMIC VISUAL RESPONSE ALGORITHM 

The flow chart for the Dynamic Visual Response Algorithm appears 
in Figure 16. The first input is simply the acceleration in units of G's 
(G = 9. 81 m/ sec [Rogers, 1978]). This input must be generated from the 
aircraft dynamics. This algorithm assumes that the acceleration is readily 
available. The acceleration is used to determine the cardiovascular 
protection variable contribution, PVC. The cardiovascular response 
model is represented by a linear transfer function which corresponds to 
the response of the pilot's blood pressure to acceleration. The calculations 
used in the transfer function are implemented with a recursion relation. 
As with all the various protection variable contributions a positive increas- 
ing value indicates an increased protection or tolerance level. 

ERIC 31 



InitiAiitation 
C » -21.4 
A > • 0.38261 
B > 1.375342 





C » -0.0947*G 
D * 0.00j634*G 
E > 0. 09831«G 
F» 0. 17453 



p 




3.32 


SEAT 


s 


0.226893 


EYE 


s 


0.087267 


Z 


» 


57.2958 



X(l) 


X 


0.0 


X(2) 


9 


0.0 


X(3) 


S 


0,0 


X(l) 


» 


0.0 


X(2) 


9 


0.0 


X/3) 




0,0 


SI 


M 


0.0 



READi X(M. DELP 



XiU • X(l) - I 



A*Y{3) » B*Y{2) ♦ C»X<3) * D«X(2) ♦ E*X{1) 






PVC » 0.0 


r— 




- 0. 0 





ERIC 



Figure 16. Dynamic Visual Response Algorithm 
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Figure 16. Dynamic Visual Response A'"orithm (continued) 
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Figure 16. Dvnamic Visual Response Algorithm (concluded) 
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The second input is the difference between the subject's actual G~suit 
pressure and the "standard curve pressure * (Rogers, 1978) (as calculated 
from the standard pressure suit valve fill schedule). This difference, 
which IS labeled DELP, is used to calculate another protection variable 
contribution, PVG. The final input is from the EMG algoi ithm and consists 
of a single protection variable contribution, PVS. This input is explained 
more completely in the EMG algorithm section. 

After ail the protection variable contributions have been calculated, 
they are summed to form the protection variable, PV, This protection 
variable is then used to calculate the blood pressure at eye level, PAE. 
The variable PAE is then used to determine the size of the visual field. 

The size of the visual field is then used to determine equivalent 
light levels for center and side CRT screens. Brightness and contrast 
signals for the center and side screens are then calculated by scaling these 
light levels. 

This algorithm assumes that all three inputs are readily accessible, 
and the protection variables are positive. This algorithm was written in 
Fortran IV and was tested on both the University of Dayton Sper ry-Univac 
70/7 Computer and the Air Force Human Resources Laboratory (AFHRL) 
Simulation Training and Research System (STARS) Xerox Sigma 5 computer 
facility at Wriglit-Patter son Air Force Base. 

3.3 FMG ALGORITHM 

This algorithm is designed to take the rectified EMG signal from 
the subject's straining muscles and generate the protection variable con- 
tribution due to straining, PVS. The variable PVS in turn will be used in 
the main program as one of the protection variable contributions* In the 
process of calculating PVS, this program also calculates the EMG variable 
PVK^ and the straining effectiveness value, K2. Both these variables, 
PVM and K2, are important because they are separate functions of time. 
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Figu*-e 17. EMG Algorithm 
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Their product determines ♦'he straining protection variable contribution, 
PVS. The < Igorithm flow chart appears in Figure 17. 

The single input of this algorithm is the rectified Electromyographic 
(EMG) signal from the subject's uody. The EMG signal is used only in 
the determina of the EMG protection variable PVM. Upon the initiation 
of a desired level of straining the counter v<iricibie flME is started. After 
TIME is initiated the straining effectiveness value is calculated according 
to the current value of TIME. 

The single output of thi • algorithm is the straining protection 
variable contribution, PVS. The variable PVS is calculated with the above 
ment ned variables K2 and PVM. Also a scaling variable, R, is used to 
keep the value of PVS between acceptable levels. Appropriate limiters 
have been written into the algorithm to keep the variable PVS positive. 

This algorithm was written in Fortran IV and was tested on both the 
Uni ersity of Dayton Sperry-Univac 70/7 Computer and the AFHRL STARS 
facility. 

3.4 PROGRAM CHECK-OUT 

Th'ise two algorithms have been combined together into a single 
program for use as software for the BIOCONAID system. The listing of 
this program, called BIO, is shown in Appendix B. 

Six different computer generated input profiles were used for 
program check-out. Three of these had combined acceleration and 
straining inputs, while three contained acceleration inputs only. The 
results of these profiles, in terms of the protection variable, PV 
(calculated blood pressure at eye level), are presented in Figures 18 to 
23. All of these results exhibit a delayed loss in PV due to acceleration 
and a delayed increase in F due to straining. In addition, these results 
agree with the mooels presented in Section 2. 
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Figure 19. PV Response to Profile #2 (No EMG Input) 
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Figure 20. PV Response to Profile #3 
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Figure 21. PV Response to Profile #4 
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Figure 23. PV Rpsponse to Profile #6 (No EMG Input) 
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6. 

SECTION 4 

STRAINING EXPERIMENTS WITH HUMAN SUBJECTS 

4. I INTRODUCTION 

The experiments and equipment development described here pertain 
to the EMC subsystem which was described earlier in the report (Section 
2.5). These experiments were done to accomplish the following goals: 

1. To evaluate possible electrode sites on the pilot's body to pick 
up straining signals. 

2. To show that a reliable, usable signal is produced. 

3. To indicate tolerance of the EMG subsystem to extraneous body 
motions. 

4. To aid in the development of electronic filters for the EMG 
subsystem. 

The experiments involved recording the EMG signals produced as the 
8'::bjccts strained in response to visual commands. The seven subjects 
were faculty members and students. All faculty members were experi- 
enced pilots. The work was accomplished at the University of Dayton and 
at the Air Force Human Resources Laboratory, Wright- Patterson Air 
Force Base, during the period of November 1978 to February 1979. 

4. 2 EXPERIMENTAL SET-UP AND PROCEDURE 

Figure 24 sho\vs the experimental setap. The equipment used and 
the experimental procedure ar^* described below. 

Electrodes and wiring 
Bio-signal amplifiers (Gain =^ 100) 
TR-20 Analog (Hybrid) Computer (EAI) to provide: 

a. Hi-Pass Filter 

b. Amplification 

c. Signal Rectifier (full wave) 

d. Lo-Pass Filter for smoothing 

e. Voltmeter - used for visual feedback of error signal 
to subject 

Strip Chart Recorder 
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Figure Z4» Experimental Setup 



Procedure; (equipment "on" and checked out earlier) 



t5 minutes) 1, 



Prepare subject s Strap-on electrodes. 
Explain M-l/L-1 straining maneuvers, 
straining will make error needle move. 



Use EMG Lotion, 
Explain that his 



(3 minutes) 2. 



Subject practices straming ind observes needle* Operator 
inserts practice commands and subject strains to bring 
needle back to zero* 



(Z minutesi V 



minutes) 4. 



First test period (Experiment No. 1 ). Electrodes on 
back. Subject relaxed and motionless. Command 
level one (3 volt) is inserted and subject strains to 
keep error zeroed. After about five seconds, 
command is removed and subject relaxes. Repeat, 
when subject ready, using commancl level two (6 volt), 

S^cont test period (Experiment No. Z). Repeat 
except subject periodically raises and lowers alter- 
nating leet {\-2 inches) during straining periods. 



minutes) 



C hange electrodes to abdomen. 



{I minutes) S. Third test period (Experiment No. 3). Repeat first. 

(Z minutes' r,. Fourth test period (Experiment No. 4). Repeat second 
except subject raises and lower alternating arms 
(^j-B inches) during strainirm period instead of 
raising fpot. 
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(5 minutes) Change electrodes to buttocks (seat), 

(2 minutes) 7. Fifth test period (Experiment No. 5). Repeat first. 

(2 minutes) 8. Sixth test period (Experiment No. 6), Repeat second 
except extraneous motion is turning head (about 45^) 
alternately each way during straining period. 

(3 minutes) 9. Remove electrodes. Equipment off. 

(33 minutes total time) 

4.3 RESULTS 

Preliminary results are presented in Figures 25-27 for three different 
subjects. This data was taken with an early technique of EMG signal pro- 
cessing. The upper trace is the commanded signal, which the subject 
attempts to match with his smoothed EMG signal (the middle trace). The 
bottom trace is the raw EMG signal. 

Although these preliminary results gave useful raw EMG signals, 
the smoothed EMG signal suffered from a low frequency drift and excessive 
noise. In addition, the smoothed EMG signal exhibited an excessive time 
lag compared to the raw EMG signal. 

Several changes in the EMG signal processing technique overcame 
these shortcomings and are discussed in Section 4.4. The results for 
experiments using five subjects conducted at the University of Dayton with 
the final signal processing technique are presented in Figures 28-33. 
These dala show a well-behaved smoothed EMG signal, free of drift or 
extraneous no.ie. 

In addition a number of runs were made with a different subject m 
the STARS simulator cockpit at AFHRL. The BIOCONAID EMG subsystem 
equipment was operated at the STARS facility to determine whether thei 
would be an excessive amount of electronic noise interference in that 
environment. These results are presented in Figures 34-40, in terms of 
the raw EMG signal and an error signal (commanded signal minus 
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smoothed EMG output). The equipment operated satisfactorily without 
shielded-cables and amplifiers. All STARS equipment, including the TV 
view screens, was operating with the exception of the cockpit instruments 
(400 Hz). 

4.4 DISCUSSION 

Early in the experiments it became apparent that a more effective 
filter was needed for the EMG signal. The smoothed EMG signal had too 
much noise associated with it. The raw EMG signal included a low 
frequency, drift, which also affected the smoothed signal. The smoothed 
EMG signal also exhibited an excessive time lag. 

Several changes were made to the EMG signal processing technique 
to improve the system. A high-pass filter was added to remove the DC 
drift from the incoming EMG signal. The low-pass filter (smoothing net- 
work) was converted to a two-stage filter. Time constants of all filters 
were adjusted for near-optimum response. The final configuration is that 
shown in I gures 13 and 24. Figures 28-40 show a smoothed EMG signal 
which IS reliable, free from drift and free from excessive noise. This 
smoothed EMG signal shows very good correlation with the straining (raw 
EMG) signal from the pilot. 

The Experimental Procedure lists Experiments 5 and 6 which were 
to be done with electrodes on the buttocks. Only Figures 25 and 26 show 
any such tests being accomplished. The reason Experiments 5 and 6 were 
discontinued is that consistently much better results were obtained with 
the lattisimus dorsi (back) muscles and with the external oblique (abdominal) 
muscles. Figure 26 shoves a comparison between Experiments 4 and 5. 
The smoothed EMG response to the 3 volt command has roughly three 
times the magnitude in Experiment 4 (abdomen) than in Experiment 5 
(buttocks). Similarly, the raw EMG traces showed a marked difference. 
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Another problem encountered with placing electrodes on the buttocks 
was the difficulty in strapping them on. This problem would be eliminated 
with the adhesive electrodob used in later tests. 

It appears the ^reat advantaj^e of seat electrodes v.oald be their 
possible use without the necessity ol mountin:; electrodes on the pilot's 
skin. Efforts to develop "seat pan electrodes should be encouraged. f 
electrodes must be mounted on the? pilot, then the abdomen or back ^^ives a 
much stronger signal. 

Experiments 1 and Z had electrodes on the back while Experiments 3 
and 4 had electrodes on the abdomen. Examination of the raw EMG traces 
m Figures Z8 through 33 indicates either site ib satisfactory but the 
abdomen usually gives a slightly stronger siL^nai. 

Experiments Z and 4 involved having the subject do extraneous body 
motions while straining. Rhythmically raising and lowering feet (Experiment 
Z) and raising and lowering arms (Experiment 4) were the specific motions. 
Experiment 6 involved rhythmically turning the head from side to side. 
There was no noticeable effect on the smoothed EMG signal due to these 
motions. However, major motions (standing up or bending iat forward) 
would saturate the signal processing equipment. Also touching or jostling 
the electrodes caused a large response. We feel that the normal motions 
of a pilot in the simulator cockpit will not adversely affect the BIOCONAID 
system. 

Two input command levels (usually 3 and 6 volts) were used to show 
that the pilot can control his level of straining with reasonable accuracy. 
Figures Z8 through 33 show the smoothed EMG signal approximating the 
command level. 

Adhesive electrodes were obtained in time to be used in part of the 
University tests on one subject (Figure 33) and on all the tests in the STARS 
facility c.t AFHRL. As might be expected, these electrodes maintained a 
more constant skin contact then the strap-on type electrode. Also, they 
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Figure 30. Subject #5 
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Expe riment # 
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h igure 3Z, Subject #7 
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Figure 35. Subject #8, STARS 
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Figure 39. Subject #8, STARS 
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are easier to mount, may be used several tunes, and the lotion in them 
may be replenished i"* dr/. Flowever, removal oi the e adhesive electrodes 
IS not pleasant since they pull hairs out. 

A difficulty with amplifier saturation in the was initially 

encountered at the STARS facility. Howe\er, v.hen the gam was reduced 
by a factor of 10, e\erythini^ operated normally. To compensate for the 
reduced uain, comoian.l su^.nal levels were .Uso reciuced to 0.3 md O.o 
volts. 

Som( very hii^h frequency (Z MMz;, low amplitude, noise pick-up was 
observed r.n an oscilloscope connected to the raw EMG signal. This may 
have been due to lae TV. Since it did not interfere with the BIOCONAID 
system operation, it Vvas concluded that electronic noise pick-up at the 
STARS facility would not be too severe. 



SECTION 5 
CONCLUSIONS ANL RECOMMENDATIONS 



5. I CONCLUSIONS 

The BIOCONAID System offers the potential to provide a very 
realistic simulation of visual dimming with straining during large accelera- 
tions. The Dynamic Visual Response Algorithm provides a simulated 
visual response in agreement with experimental observations and known 
physiology* The EMG Algorithm provides an input to the visual response 
which is also in agreement with experimental observations of straining 
during large accelerations. 

The EMG subsystem provides a smoothed, reliable straining signal, 
not adversely affected by extraneous limb movements. Either abdomen 
(external oblique) or the back (latissimwas dorsi) muscles are suitable as 
electrode sites. 

All parts of the BIOCCNAID System have been carefully checked-out 
with the exception of the Visual Field Model. Equipment limitations at the 
STARS facility prevented the full implementation of the Visual Field Model 
during the contract period. Thus final determination of system effective- 
ness must await complete system implementation and evaluation in an 
appropriate simulator. A test plan for such implementation and evaluation 
IS contained in Appendix A. 

5 . 2 RECOMMENDA TION 

It is strongly recommended that the BIOCONAID System be imple- 
mented, demonstrated and evaluated in an appropriate Air Force simulator 
with experienced pilots. Although the system appears to provide a realistic 
simulation of visual dimming vhile straining, human pilot evaluations are 
necessary to determine system effectiveness. It is expected that such 



iniplenutntatLon would invol'vc tnie acijustnient ol bon.e algorithm i>ains, 
thoreby enhancing the by stem. 

It lb cilbc^ recummended that sonx* medical research be conducted into 
the phybiolot^ical etfectb of sustained straining m a 1-G environment. It 
W'.ulfi bo extrensely ubetul to identify the level of straining which causes 
bignificant detrinu-ntal cardiovascular effects, as well as the level of strain- 
ing which causes negligible detrimental effects. This medical research 
should bemn with a literature review of experimental and historical data. 
If this review floes not yield definitive aiiswers on the above straining levels, 
n»*vv experin^ents on human subjects or animals should be conducted to yield 
th»' dt sired results. This infornmtion is necessary for the safe implementa- 
tion "I the lUOCCNAID System in a training environment. 
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APPENDIX A • - 

TEST PLAN 
BIOCONAID 

(Bionic Control of Acceleration Induced Dimmirv^) 
SYSTEM TEST AND EVALUATION 

A. 1 INTkODLC nON 

This tcfet.plan pertains to the perf ormance'evaluation of the final 
design of a system for bionic control of acceleration induced diuimmg 
(BIOCONAID^ Present advanced flight simulators provide a system*for 
dimming the pilot's visual scene to simulate partial blackout when the 
aircraft normal accd'ieration exceeds an appropriate level. The new 
system goes one step beyond and provides for control of the dimming 
rough the pilot's muscular straining action (M- 1 or L~ 1 maneuver). 
The BIOCONAID system has been designed by the University of Dayton 
under Contract F336 1 5- 77-- C- 0080 with the Air Force Human Resources 
Laboratory. This test plan has been prepared to fulfill a contractua- 
requirement. 

The plan involves installing the BIOCONAID system in the SAAC 
F-4 Simulator at Luke AFB, conducting the tests using experienced 
Air Force pilots, removing the BIOCONAID system and writing a report. 
The intention is to minimize interference with normal training use of the 
simulator during the evaluation. A seven week period is provided for 
installation, interfacing and check-out followed by a one-week test 
period. Disassembly and removal will be accomplished in one day. 

The specific time scl>cdule for the tests is dependent upon the date 
Air Force go-y4head is given and upon other Air Force instructions. The 
University of Dayton anticipates this evaluation will be conducted under 
an additional contractual arrangement. The following section contains 
the test objectives* r 
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A- 2 TEST OBJECTIVES 

A.Z.I Objective A 

To evaluate the performance of the final design , specifically, 
IS the visual scene dimming reliably and realistically controlled by pilot 
straining'' To answer this question, the system must be observed during 
ion^i-term, hiuh-acceleration turns. In addition, however, observations 
must be made during rapidly changing mariffuvering accelerations as well 
as during periods gf low dynamic acceleration. The'se latter observations 
will be for the purpose of determining whether the final design is free from 
inappropriate respcnses. All tests will be done with and without the electro- 
myographic feedback to allow comparisons to be made. 
% 

A.Z.2 Objective B 

To conduct other tests on the final system . These "other 
tests" include further evaluation of several EMG electrode types and elec- 
trode sites and further evaluation of system tolerance to extraneous motion. 
These tests can be done with minimal increase in test time beyond the 
Objective Attests and they should be very significant since they will involve 
Air Force pilots in actual simulated flight. Pilot tolerance for the electrodes 
is the most important question to be answered. Electrodes to be included 
in these tests are: 

• strap-on, reusable many times 

• adhesive, disposable 

• electrode belt 

A, 

• seat pan electrodes (if cwaH^tble) 
Electrode sites included are: 

• abdomen (external oblique) 

• back (latissimus dorsi) 

• buttL k (with seat pan electrodes) 
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A. i TEST PROCEDURE 
A» 3» I General 

During a ^ne-week test period a numb^ir ol experienced Air 
Force pilots will act as test subjects to perform simulator iiights aimed at 
accomplishing the test objectives. xae planned number participating is 
ten (10), although more would be welcome. Each pilot will be briefed, 
electrodes will be installed and a short practice period with the BIOCON- 
AID system will be provided. Next, the actual test maneuvers will be 
accomplished. These maneuvers consist of long-term, high accelera- 
tion turns in bo^h directions, short-term, rapidly changing acceleration 
^arns or aerobatics and normal operations near one-G, such as landing 
pattern. Details of the test protocol are given below. Finally, after the 
flight, the electrodes will be removed and the pilot will be debriefed. 
Each pilot will provide a brief written observation report which will be 
included in the final test report. 

rhe actual time -in-simulator for each pilot is estimated as 
30-45 minutes and the total time per pilot is estimated as 1—1-1/2 hours. 
The time of day for the tests will be such as to minimize interfering with 
^tbe normal simulator training schedule. Once the BIOCONAID equipment 
1b thstalled, it can be used or not used as desired. Merely placing a switch 
in the OFF position returns the simulator to its original configuration. 

A. 3. 2 Specific Test Protocol ^ 

Time (min)/Cum . Activity 
13/13 Pilot Briefing: 

- Review blackout phe*>^.xiena 

- Review M- 1 / L- 1 maneuver 

- Explain BIOCONAID system 

- Explain test objectives 

- Describe test procedure 
• The need for electrodes 
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Install electrodes according to the following table: 



No, Pilots Type Location 

1 Strap'On back 

1 Strap- on abdomen 

2 Belt waist 

3 Paste-on back 

3 Paste - on abdomen 

If seat pan electrodes are available, tliey will be 
used with all pilots (i. e. for part of each pilot-test 
period ). 

l(y/25 Pilot enters and starts simulator, climbs to 5000 

ft altitude and practices turns to familiarize him- 
self with BIOCONAID. 

20/45 Test Period. Pilot does the following maneuvers 

at 5000 ft, Mach . 8 on instruction from controller 
(roll level after each turn): 

- A 360° turn right with 60° bank rngle (ZG) 
' Repeat to left without EMG feedback 

^ A 360° turn right with 70° bank angle (2. 92 Gj 
twist body and head to look over right shoulder 
Repeat to left without EMG feedback, twist 
body and head to look ovej- left shoulder 

^ A 360° turn right with 80 bank angle (5. 76G) 

- Repeat to left without EMG feedback 

- A series of alternating right and left turns of 
45° with bank angle 70°, rolling directly from 
one turn into the next. Eight turns will be 
done flast four without EMG). 

A Chandelle right followed by a Chandelle 
left (no EMG) 

- A Lazy- Eight (no EMG when turning left) 

- A rectangular landing pattern 

5/rQ Pilot shuts down simulator, c^eplanes, and electrodes 

are removed. 
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Pilot is debriefed and gives written observations: 



- Tolerance for electrodes. 

- Was energy expenditure to avoid blackout 
realistic '> 

- Compare with and without EMG as regards 
realisnn of simulated flight. 

- Was the onset of blackout and the recovery- 
through straining realistic 

- Was the beneficial effect of the G-suit realis- 
tically simulated? 

- Tolerance of system to bodv twisfing. 

- Further comments or recommendations. 



Figure 41 shows a block diagram of the integrated system 



including pilot, aircraft, G-suit, cardiovascular model, straining model, 
and visual scene model. The numbers and names used below for recorded 
variables refer to Figure 41. 

Both digital and analog signals v/ill be recorded. It is 
anticipated that digital variables will be transferred to tape and printed off-- 
line. The analog recorder will be turned on and off by one of the operators 
as the test progresses. 

Teletype input to the simulator computer will result in 
appropriate printed idertification of data pertaining to each maneuver. 
When a new maJieuver is due, the teletype input is changed by the operator. 



A»3.3 Recordings 



Table A-1 lists variables to be recorded together with type 
(digital or analog) and frequency of observation. 



A. 4 TIME SCHEDULE 



(Note: T = First test day, designated by AF) 



TIME 



EVENT 



T-49 



UD personnel (3) visit Luke AFB for test planning 



T-39 



UD ship or transport equipment to Luke AFB 
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TABLE A-1 



TABLE OF RECORDINGS 

Each Each 
Type Once Maneuver Second Continuous 



Date, hour, subject D 

Electrode Location D 

Electrode Type D 

Maneuver D 

Time D 

Roll Angle (^) D 

Heading Angle (t) D 

G-level (N^) D 

Raw Straining Signal A 

Processed Strain Signal A 
(EMC) 

G-Suit Protection Value D 

Cardiovascular Model D 
Output, PV^ 

Straining Protection D 
Value, PVg 

Net Protection Value, PV D 

Dimming Signal (Output A 
of Visual Scene Model) 

Visual Scene Light A 
Intensity 



X 
X 
X 
X 



X 



X 



X 
X 
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G(t). 



SEAT ANGIE 



tARDIOVASCULAR 
MODEL 



G(tl 



G SUIT MODEL 



AIRCRAFT 
DYNAMICS 




VISUAL FIELD 
MODEL 



STRAINING 
MODEL 



PROJECTI 
SYSTEM 



VIEW SCREEN 

\ 

INSTRUMENT 
PANEL 




STRAININI 
SIGNAL 



CONTROL COMMANDS 



-H- 



Figure A— I System Integration Diagram 
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T-37 UD engineer/programmer arrives Luke AFB, began interfacing 
equipment and software 

T-27 Interfacing complete 

T-17 Piloted system checkout starts (one AF pilot needed) 
T- 10 BIOCONAID system check-out complete 

T-5 Forijial demonstration of systenn and AF approval to proceed 

with test 

T Testing starts 

T-f6 Testing ends 

T+7 Equipment removed, simulator restored to original configuration 

T+67 Draft copy of final test report to AFHRC 

A. 5 IMPACT ON SAAC SIMULATOR SCHEDULE 

It is expected that the simulator designated for these tests will be 
already fully utilized during normal duty hours for training activities. 
Consequently, the training usage should retain first priority and the test 
program activities would proceed on a *'time available'^ basis. Hope- 
fully, evening and weekend hours would be available regularly during 
the 44 day set-up and test period. It is estimated that ?.n average of 
two hours operating time per day would be required. ^ 

Care will be taken to insure that installation of the test equipment 
and software is done in such a way that it does not affect the normal 
training use of the simulator. This means that simple switches and 
teletype entries will disarm all test provisions and leave the simulator 
in its original operating mode. 

Also, the test provisions will be installed so that they may be 
rapidly removed leaving the simulator in its original configuration. 
Complete removal of test provisions from the site wiT require less 
than one day. 
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A. 6 AIR FORCE SUPPORT REQUIREMENTS 



A. 6. 1 Equipment 

SAAC (F-4) Simulator (8^' operating hours) 

Includes: - digital computer capable of handling software 

- one analog/digital p .rt 

- light dimming capability 

Four channel strip chart recorder. 
A. 6. 2 .A per sonnel 

1 Project Officer - liaison with base (part-time) 

1 Simulator operator /maintenance technician - to operate 
^ ^ ^ simulator (if required), monitor test installation (part> 

time), answer questions 

10 (or more) voluntary pilot subjects (1-1/2 hrs each) 

1 pilot for "piloted system check-out'' prior to tests (4 hrs) 

A. 6. 3 Work Space 

Combination off ice /wo rkroorn near simulator. 
Access to simulator. 

A. 6. 4 Security 

Secure area for UD equipment (table top analog computer, 
meters, oscilloscope, oscill?'«'r, bio-amplifiers, cables, 
tools). 

Security passes for four UD project workers. 
Car pass for rental car. 

A. 7 UNIVERSITY OF DAYTON SUPPORT 

A. 7. 1 Equipment and S'lpplies 

Oscillator 
Bio-Amplifiers (6) 
Cables and wiring 
EMG electrodes 

/ if ^ » 
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A- 7. 2 Personnel ^ 

Z^engineers (one visit per week) 

2 technicians or students (full-time on ^^^.m 

A*J* 3 Transportation, Housing and FoO':^ 

The University will furnish transportation for its personnel 
and equipment ai^y^'ell as housing and meals for its personnel. 

A. 8 REPORT^REQUIREMENTS 

/ 

A report] will be prepared as prescribed by Air Force instructions- 
A draft copy will be furnished ,o AFHRL oO days after completion of the 
tests. The report will include the recorded data, the data analysis, the 
pilot comments, conclusions and recommendations, and a discussion. 




83 



APPENDIX B 
ALGORITHM LISTINGS 



00010 L 

oooao 

000 30 
'\ 015134^1 
60050 
00060 
00070 
00080 C 
00090 C 
00100 C 
00110 c 
00120 C 
00130 C 
00140 C 
00150 C 

00160 •: 

^ 00170 C 
, 00180 C 
00190 C 
00200 C 
00210 C 
00220 C 
00230 C 
00240 C 
00250 C 
00260 C 
00270 C 
00280 C 
00290 C 
00300 C 
00310 C 
00320 C 
00330 C 
00340 C 
00350 C 
00360 C 
00370 C 
00330 C 
00390 
00400 
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TFDPSPC.BID IiRVID QUflM 

rCMMDN'INPUTS/X<6' 'DEI P.PV^ 

CDMilDN/DUTPUT.^BPS ' BPC ' C T : . C TC 

DIMENSION Y'3;' 

PERL N.NFM.K2.LLS.LLC 

CB = 1.0 

CC = 1.0 



INITIRLIZRTION 



G = i^fllN DF CRPDIDVRSCULRP TPRNSFEP FUNGTIDN. 
ft.B»C,D'E = TIICTIN CD-EFFICIENTS DF TPf^NSFEP FUNCTION. 
F ^ OFFSET mSTRNCE OF THE CPTIC DISC FPOM THE CENTPRL 

VISURL RXISdO DEGREES FOP PIGHT EYE> 
P = OFFSET DISTRNCE OF THE OPTIC DISC FROM THE CENTPRL 

VISURL RXISaiO DEGREES FOR THE LEFT EVE) 
SERT = ^ERT RNGLE CSTRNDRPD IS 13 DEGREES^ 
EYE = OFFSET RNGLE OF EYE FROM PHYSIOLOGIC Z-RXIS 

CTYPICRLLY* EYE = 5 DEGREES^ 
Z = -CONVERSION FRCTOP FDR DEGREES TO RRDIRNS CDNVEPSIDN. 
Xfhl) = G => DIMHNSIDNLESS MULTIPLE OF GPRVITY CONSTRNT. 
•i<i> = OUTPUT OF CRRDIOVRSCULRR TRRNSFER FUNCTION. 
DELP = DIFFERENCE BETWEEN STRNDflPD PRESSURE SUIT VRLVE 

FILL SCHEDULE RND THE RCTURL SUIT PRESSURE. 
PV = TOTRL PROTECTION VRLUE. 
PVC = CRRDIOVRSCULRR PROTECTION VRLUE. 
PVS = ELECTROMYOGRRPHIC PROTECTION VRLUE. 
PVG = G-SUIT PROTECTION VRLUE. 

PRE = CRLCULRTED BLOOD PRESSURE RT EYE LEVEL. 

SI = POLRR POSITION RNGLE IN VISURLFIELD MERSURED 

COUNTERCLOCKWISE RBOUT VISURL RXIS FROM RIGHT 
HORIZDNTRL RXIS. < 0, LE. SI . LT. 2»PI RRDIRNS.;- 

INITIRLIZRTION FOR DYNRMTC VISURL RESPONSE ALGORITHM 



G = -21.4 

fl = -0.882609 

Inn 

■* V/ ^ 
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0 0-^ 1 0 
00420 
00430 
0044 0 
00450 
0046 0 
00470 
00480 



P 1.37^342 

L = 0'?47 ♦G 

E = 0. U*?&51 ♦b 

F = 0.17453 

P = 3. 3d 
ZEHT = 0, •326893 

E'r'E - 0. 0372665 



0049 0 




2958 


0050U 


X»l' = 


0. 0 


0051 0 


/; ^ 2 ' = 


0 . IJ 


00520 




0. 0 


0 05 3 0 


, 1 . = 


»"'. u 
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0 065 0 
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00680 
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00700 
00710 i 
00720 f 
00730 
00740 ( 
00750 ( 
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00780 C 
00790 L 
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00810 C 
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00830 C 
0084 0 C 
IJ0050 L 
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00890 
00900 
00910 
009r0 



INITIRLIZhillDM FOP EMG HLGDPITHM 

H = 0 

HNM = 0 

NFH = 5 

P = 1.5 



1 PEftD^99* 101 -'X* 1 ' 'EMb'DELP 

101 FDPMhTi 3'F14.5' ■ 
/>1 • = ;>1 » - 1 



NNM ^- NNM ^ 1 
EMG HLGDPITHM 

DETECT :TPhINING 

IF^EMG.GT.0.7> GD'TD 15 

HD ITPRiNING SD LET PVS = 0*0 

H = 0 
PVS = 0. 0 
GD TD 21 

ITPftlNING DETECTED == ITfiPT CDUMTEP 
CRLCULRTE RND LIMIT PVM 



0^^30 C 
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X V/ u 



00940 

00950 
00960 
00970 
00980 
00990 
01000 
01010 
01020 
01030 
01 040 
01 050 
01060 
01070 
01080 
01090 
01 100 

OHIO 

01 120 
01130 
01 140 
01150 
01 16G 
01170 
01180 
01 190 
01200 
01210 
01220 
012 0 
01240 
01250 
01260 
01270 
01280 
01290 
01300 
01310 
01320 
01330 
01340 
01350 
01360 
01370 
01380 
01390 
01400 
01410 
01420 
01430 
01440 
01450 



15 

C 
C 
C 

C 



10 
11 

12 

C 
C 

20 



21 

C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 



c 
c 
c 
c 
c 
c 
c 
c 
c 



PVH = <22.39*EI16)-15.67 
IF<PVM.GT.51.5^ PVM = 51.5 

DETERMIME EFFECTIVENESS DF STPfilNING DUE TD TIME 
FftCTDP <K2) 

TIME =M.-^NFM 

IFaiME.LE.2. 0^ 3D TD 10 
IF<TIME.LT.4.0) bD TD 11 
IF<^TIME.LT.5. 0> 6D TD 12 
K2 = <16. 0 - TIME^ /I 1. 0 
GD TD £0 

K2 = 0.5 ;' 

GD TD 20 

k2=.£5*TIME 

GD TD 20 

K2 = 1.0 

DETERMINE FHSD LIMIT PVS 

PVS = K2 ♦ PVM ♦ P 

IF<^PVS.LT. 0. O^PVS = 0. 0 

N ^ H 1.0 

CDNTINUE 

♦♦♦♦ » ♦»»»»»»♦♦♦♦♦*»»»»»»♦»»»»»»»♦♦♦■»♦♦♦♦»♦»»♦♦♦♦♦♦♦♦♦ 

EMD DF EMG BLGDPITHM 
♦♦♦♦♦♦♦»»♦»♦♦♦♦♦»»»»»»»♦♦♦♦»»♦♦»»»»»»♦»»♦♦♦-»♦♦♦♦♦♦♦♦ 

»»»»»♦»♦♦»»♦»♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦»♦»»♦»♦♦♦»»■»■♦»»♦♦»♦♦♦ 

CftPDIDVRSCULftP RESPDNSE RLGDPITHM 
♦♦♦♦♦♦♦♦-»♦♦♦♦♦♦♦ » »» » »♦♦♦♦♦»«■»»♦»»»»»»»»»»♦»♦♦♦»»♦♦»♦♦ 

CftLCULRTE PVC FPDM CftPDIDVRSCULRP TPfiNSFEP FUNCTION 
»♦♦♦♦♦♦♦♦♦♦♦♦♦♦»♦♦♦♦♦♦♦♦♦♦ »»»»»»♦»»»♦» »♦♦♦»♦♦♦♦»♦♦♦♦♦ 

Y^n = R*Y< 3'* + B*Y<:2> + C*X< 3) + D*X (.a;" + E*X<0 

PVC = 120 + Ya> 



LIMIT PVC ==> NDT TD BE LESS THRN 0 
IF< PVC. LT. 0.0 ^ PVC = 0.0 
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G^SUIT MODEL 

CflLCULflTE PVG ==>NDT TD EXCEED 42.8 
PVG=0. 

IFcDELP«GT.-80> PVG = <42.8^(DELP 30> /80>^X<.1> 
IF^PVG.GT.42.8; PVG = 42.8 

♦^♦♦♦^^♦♦♦♦♦♦♦♦♦♦♦♦^^♦♦♦»»»»»»»»»»»»»»»»»»»»»»»»»»»#^ 

CFH^CULfiTE PV ~ TDTflL PRQTECTIDN VfiLUE 

»»»»»»»»»»»»»» 

PV = PVC PVG PVS 

♦♦♦♦♦♦♦♦♦♦♦»♦»♦»♦♦♦♦♦♦♦♦»♦♦»♦»»♦♦»»»»»»»»»»»»»»»»»»»» 

CftLCULfiTE PRE 

♦♦♦♦♦♦♦♦♦♦♦»♦♦♦♦♦♦♦♦♦'»'♦♦♦♦♦♦♦»»»»♦♦♦»»»»»»♦♦»»»»»♦♦♦- 

PRE = PV - <0.77*3i:i.48*CDS<.SEftT ~ EYE>*X<;i>> 

IF PRE IS LESS THBN 37.3, VISION CDLLfiPSES 

ir <:PV.GT.37.3> GD TD 7 
LLS = Ci. 0 
'^LC =0.0 
WPITE<1,88> 

WRITE ( 98 f 87^' NNM^ X < 6," EMG' PV 
FDPI1fiT'5X' "VISION CDLLflPSE") 
GD TD 32 

♦♦♦♦»♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦»»»»»♦» »♦♦♦♦♦♦♦♦»♦♦♦♦♦♦♦♦♦♦♦• 

CfiLCULfiTE SENSED VISUfiL fiNGLE (RRDEG, RLDEG> 
OlE CftLCULfiTE THE VISURL L IMI T, RRDEG»fiT fiNY 
fiNGLE 'SI, USING ONLY THE RIGHT EYE. 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦^iV^ 

Ul=l .91986*'?V-37> / <PV-£0) 
FI = 3.14159 

IF 'SI.Ll .Pl^a. 0,> SY = :i 
IF CI.GE.PI '2. Ci> SY = PI-SI 
TF <SI.GE,3.*PI/2. 0> SY = -2.*PI-»-SI 
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01980 C 

01990 C HDU -PI '2. .LE. :Y.LE.-^P1 '5. 
02000 C 

02010 ftH=F^i:n'; • . i745 3-^:v » 

02020 EB=-:OPT^- OSI^'CQ: ' . ir453-»^:V' '♦♦2-- 0::l-»^»'i^^2 • 

02030 PIP- = FIH PB ' 
02 040 r 

02 05 0 C ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

02060 C 

02070 c r Man IE coppect :en:ed vicuhl hngle hmd corvv'EPT 

02 030 c TD degpee: 

02090 L 

02 1 0 0 C ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦^♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦^^^^^^^^^ 

02110 L 

02120 PPDEG = Z ♦ PIP 

02130 LJ-.: = • P^'-l 1 U- • 4. 0. 

02140 LLC = »PV-3S. • 5.0-»^10. 

02150 IF^LLI-GE- 10- • LLI = 10. 

02160 IF^LLC-bE. 10. • LLC = 10. 

02170 IF^LLZ.LE. 0. ' LLI = 0- 

02180 IF^LLC.LE. 0. • LLL=0. 

02190 32 bp:=ll:kb 

02200 bpc = llc ^ cb 

02210 CTT = LL: ♦ CC 

02220 - CTC = LLC ♦ CC 

02230 IF'LLC-LE. 0. > LLC = 0. 

0224 0 L 

02250 C ♦♦♦♦♦♦♦♦♦♦^♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦^ 

02260 C 

02270 C OUTPUT lENIED VIZUftL HiibLE HND UPDhTE RPPRVS 
02280 C 

0229 0 C ♦-^♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦^ ♦♦♦♦♦♦♦♦♦♦♦^♦♦♦♦♦♦^ 

02300 C 

02310 3 X'6> = X» 1 » ^ 1 

02320 PPIhfT X^6' •EnG»PV»LL-»LLC 

02330 WRITE <98»87>mn»X '6* •EflGf PV 

02340 37 FDRMRT M2!t 3Fr». 3 ' 

02350 33 X O • = X » 2 ' 

02360 X f 2 • = X • 1 • 

02370 XM' = 0.0 

02580 Yf3» = V'2> 

02390 Y't* = V • 1 • 

02400 V a ^ = 0- 0 

02410 C / 

02420 C RUN FINIiHED ' 
02430 C 

02440 IF<mn.NE.99^GD TD 1 

02450 :tdp 

02460 END 
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APPENDIX C 
THE M-1 AND 1 STRAINING MANEUVERS 

C. I M-1 MANEUVER 

''Pilots commonly refer to the M-1 Maneuver as the 'grunt' maneuver 
since it approximates the pnysical effort required to lift a heavy weight. 
The M-1 maneuver consists of pulling the head down between the'^shoulder s, 
slowly and forcefully exhaling through a partially closed glottis , and simul- 
taneously tensing all skeletal muscles. Pulling the hea'^. downward gives 
some degree of postural protection (shoitens the vertical head-heart 
distance): intrathoracic pressure is increased by strong muscular expiratory 
efforts against a partially :losed glottis; and the contraction of abdominal 
and peripheral muscles raises the diaphragm and externally compresses 
capacitance vessels. Tor long-duration G exposures, the maneuver must 
be repeated every 4 or 5 seconds. When properly executed, the exhalation 
phase of the M-1 results in an intrathoracic pressure of 50 to 100 mm Hg 
which raises the arterial blood pressure at head level and 'hereby increases 

toler?nce at least 1.5 G. The inspiratory phase of the M-1 maneuver 
must (be done as rapidly as possible, sic. ). " 

C.2 L-J M/aNEUVER 

"The L- 1 maneuver is similar to the M-1 maneuver except the air- 
crew member forcefully attempts to exhale against a co mpl etely closed 
glottis while tensing all skeletal muscles. Using either maneuver the pilot 
obtains equal protection, i.e. 1. 5 C/ greater than relaxed blackout level with 
or without the anti-G suit. In a 1*^72 study, subjects Vvearinq anti-C/ suits 
and performmu either the M-1 or L-1 straininu maneu\er were able to 
maintain ad -equate vision durinu centrifuue exposure of 4 9 C/ for 45 sfMrtrnd-,. 
Higher and lonu<^r runs ha^e not been attempted. Howe\'or, it is important 
to note (a word of caution) that lorcetully exhalinu against a closed glottis 




1 



v M'thout vii^orous skeletal muscular tensing (Valsalva maneuver) can reduce 
t-C tolerance and lead to an episode of unconsciousness at extremely lov^ 
G leveU. Therefore, instruction and training on the proper method of 
performing th<3se straining maneuvers is essential." 

(Source: Gillingham, K. K. emd Krutz, USAFSAM AR 10-74, Brooks Air 
Force Base, Texao. ) 
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APPENDIX D 

THE SOLUTION ALGORITHM FOR THE CARDIOVASCULAR MODEL 



A digital algorithm is used to represent the cardiovascular transfer 
function. 

21.4 cos 0 (1 + 5. 31s) 

^ (D. 1) 



1 + 3. 23s + 5. 17 



s 



This algoiithm employs Tustin's m^ethod^ for computational efficiency. 
For a given recursion period T, determined by the simulation computer 
frame rate, the output y(k) given by the following equation. 

- 6.46T + 20. 68 2T^ - 41.36 
y(k) = -y(k-2) — - y(k-l) ~ 

T +6.46T + 20. 68 t" ^ 6. 46T + 20. 68 

2 2 
„ , + T - 10.6 2T „ ,^ 2T 
+ x(k 2) — 2 +x(k-l)— r 

T -'- 6.46T + 20. 68 T +6.46T + 20.68 

T^ + 10. 26T ^ 
+ x(k)~ (D. 2) 

T + 6.46T + 20. 68 

Here, y(k) represents the output PV^ at time tj^, y(k-l) is PV^ ''k 1 ~ *Tc'^ 
and x(k) is G(t ). A complete listing of coefficients for different values of 
T is given in Table D-1. The coefficients are derived with a transfer 
function gain equaLto 1» 

. A hypothetical aircraft maneuvering proiiie is shown in Figure D-1. 
The predicted eye level protection value (PV) due to the cardiovascular 
response model is shown in Figure D-2, The PV values are derived using 
T = » 02 and the hypothetical profile as a driving function. 
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TABLE D-1 
Listmg of Coefficients for P^■ ^ Simulation 



T 


A 


B 


C 


D 


E 


0. 02 


0. 9876 


-1.9875 


-0.0102 


0. 0000 


U . U 1 U ^ 


0. 04 


0. 9753 


-1.9750 


-0. 0202 


0 . 0002 


U . U Z 1>4 


0.06 


0.9632 


- 1. 9625 


-0. 0301 


0. 0003 


U. U3i)4 


0. 08 


0.9513 


-1.9501 


-0. 0398 


0. 0006 


U* U4U4 


0. 10 


0.9394 


-1.9376 


-0. 0493 


0. 0009 


n n c r\ '? 
0 . ObuZ 


0. 12 


0.9278 


-1.9251 


-0. 0587 


0.0013 


U. ObUU 


0. 14 


0. 91t3 


-1.9127 


-0. 0679 


0.0018 


0. 069 ' 


0. Ifa 


0.9049 


-1.9002 


-0. 0770 


0. 00Z4 


0. 079 3 


0. 18 


0. 8937 


-1.8878 


-0. 0859 


0. 0030 


0. 0889 


0. 20 


0. 8826 


-1.8753 


-0.0947 


0. 0036 


0. 0983 


0. 22 


0. 8717 


-1. 8629 


-0. 1033 


0. 0044 


0. 1077 


0. 24 


0. 8609 


-1. 8505 


r\ 1 1 1 O 

- 0 . 1118 


U . Wdl 


n 1 1 A Q 

U . 1 1 D 7 


0. 26 


0. 8502 


-1, 8382 


-0. 1201 


0. 0060 


0. 126 1 


0. 28 


0. 8397 


-1. 8258 


-0. 1283 


0. 0069 


0. :352 


0. 30 


0. 8293 


-1.8135 


-0. 1363 


0. 0079 


0. 1443 


0, 32 


0.8191 


-1.8011 


-0. 1442 


0.0090 


0. 1532 


0. 34 


0. 8089 


-1.7888 


-0. 1520 


0.0101 


0. 1621 


0. 36 


0. 7990 


-1.7766 


-0. 1597 


0.01 12 


0. 1709 


0. 38 


0.7891 


-1. 7643 


-0. 1672 


0.0124 


0. 1796 


0.40 


0.7794 


-1.7521 


-0. 1745 


0. 0137 


0. 1882 


0. 


0. 7698 


-1. 7398 


-0. 1818 


0. 0150 


0. 1967 




0. 7603 


-1. 7276 


-0. 1889 


0. 0163 


0. 2052 


O.^o 


0. 7509 


-1.7155 


-0. 1959 


0.0177 


0. 2136 


0.48 


0.7417 


-1. 7033 


-0. 2027 


0.0192 


0. 2219 


0. 50 


0. 7326 


-1.6912 


-0. 2094 


0.0207 


0. 2301 



y(k) = -y(k-2)A - y(k-l)B + x(k-2)C +f:(k-l)D 4 x(k)E 
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Figure D-1 Hypothetical Aircraft Maneuvering Profile 
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